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A new approach to model the biomatter dynamics based on the field theory is presented. 
It is shown that some well known tools in field theory can be utilized to describe the 
physical phenomena in life matters, in particular at elementary biomatters like DNA 
and proteins. In this approach, the biomatter dynamics are represented as results of 
interactions among its elementary matters in the form of lagrangian. Starting from the 
lagrangian would provide stronger underlying theoretical consideration for further ex- 
tension. Moreover, it also enables us to acquire rich physical observables using statistical 
mechanics instead of relying on the space-time dynamics from certain equation of mo- 
tions which is not solvable due to its nonlincaritics. Few examples from previous results 
are given and explained briefly. 
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1. Introduction 

Presently, the human understanding of living things are still very far from perfect 
when compared to the elementary particle physics. In elementary particle physics 
all elementary constituents of matters and its interactions are relatively well known. 
They have uniquely and successfully been modeled in certain mathematical repre- 
sentations in a unified manner. Though the model, namely the standard model of 
particle physics, is still lacking from some major issues like the massive neutrinos, 
the origin of symmetry breaking leads to mass differences between particles. How- 
ever those problems do not reduce its success story of modeling the huge numbers 
of phenomena in our nature in a single framework. This success is merely based on 
the principle of symmetry and a great tool called as field theory The symmetry 
has mainly been introduced to restrict the model under the basic assumption that 
the theory should be invariant under time and spatial changes. On the other hand, 
field theory plays an important role across the concept, from the initial lagrangian 
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till the calculation of real phenomenological observations ^El. 

Another notable reason of the success story in particle physics is the exploration 
of elementary matters. The particle physicists have for decades investigated the 
most fundamental building blocks of our nature. This effort ends up, at least for 
time being, at 16 elementary particles with the additional one hypothetic Higgs 
particles, currently investigated at the Large Hadron Collider (LHC) experiment at 
CERNlal It finally can be concluded that, constructing a unified model is getting 
easier if one considers only few numbers of a unified framework is possible when one 
goes to the elementary scale with relatively small number of interacting matters. 
Even in the scale of nuclear matters at ~ 10~ 15 m which involved hundreds of 
nuclear species, there is no reliable and unified model to describe its phenomena. 

In contrast with the 'phenomenal progress' in particle physics, living creatures 
involve much more complex systems. Although all matters should be formed from 
the above 16 elementary particles, the living characteristics emerges due to the 'el- 
ementary' living building blocks with much larger sizes. If one treats the protein, 
DNA / RNA, etc as the elementary biomatters, then we are faced with the inter- 
actions of 'point particles' with the square-rooted size of quarks in particle physics. 
Note that the quark size is ~ 10" 18 m, while protein or DNA is at nanometer scale, 
i.e. one order larger than an atom size. 

Now let us consider the protein or DNA as the elementary biomatters. Although 
its number is still huge compared to the 16 elementary particles, they are the most 
fundamental building blocks of living things have been observed. The problem is 
then turned out to the applicability of particle physics methods. The main issue is 
the scale difference, in particular for using field theory. Because in the field theory 
one has a basic assumption that each interacting particle is treated as a point 
particle without any geometrical shapes. This is clearly not the case for biomatters. 
However, in physics one always considers the relative size rather than the absolute 
ones. If one can make use of field theory on the particle size matter to explain the 
observed nuclear phenomenology (all elementary particle observations are through 
the nuclear, i.e. hadronic, processes), then it is also plausible to apply it on the 
biomatters to describe more macroscopic ensemble of living things. Alternatively, 
one can also bring an ad-hoc assumption to model the macroscopic ones using much 
more microscopic, but existing modelization. The second approach is commonly 
adopted among nonlinear theorists to deal with nonlinear phenomena like sea wave 
using nonlinear Schrodinger or even relativistic Klein-Gordon equations. 

Assuming that the biomatter interactions can be modeled in certain lagrangians, 
one can further derive not only the governing equation of motion (EOM) describing 
the mechanical behavior, but also extract relevant physical observables along the 
time or temperature evolution through statistical mechanics methods It should 
be remarked that the method is reliable only if the system is at an equilibrium state. 
Similar to the case of particle physics, it means the obtained physical observables 
are valid only for a particular period of equilibrium time. 
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In this paper, few examples of modeling the biomatters using field theory ap- 
proach are presented. Beside discussing the derived EOMs and its behavior, some 
interesting results from statistical mechanics calculations are also given. 

The paper is organized as follows. After this introduction, two examples of con- 
structing lagrangians for the case of DNA and protein dynamics are presented. It 
is followed by the discussion on the statistical mechanics calculation to extract ap- 
propriate physical observables is given in Sec. 3. Finally the paper is ended with a 
short conclusion. 

2. The dynamics of biomatter 

The dynamics of biomatter has been worked out since long time ago. In most variant 
of models, the dynamics are described by certain EOM. For protein the pioneering 
works were conducted by Davydov et.al. , while Yakushevich, Peyrard-Dauxois- 
Bishop et. al. applied it into the case of DNA system l ^ l ^ l ^ l M os t; of works have 
deployed the classical Hookes force of pendulum to construct the nonlinear EOM 
with a solitonic solution of oscillator harmonics. 

As already mentioned in the preceding introduction, this approach obviously 
has disadvantages of lacking of first principle which subsequently leads to no clue 
for expanding the theory to incorporate more complicated dynamics. 

There are some improved modelings yet based on the oscillator harmonics sys- 
tem, but in the form of hamiltonian. The new interactions representing either inter- 
nal or environmental effects are incorporated as potential terms. Using the hamil- 
tonian, one can derive analytically the physical observables through its partition 
function using statistical mechanics. These methods predict interesting phenomena, 
for instance the anomaly of negative heat capacity in the case of monomer with an- 
harmonic oscillation effects in thermal bath. The anomaly has been observed at 
certain region of lower finite physiological temperature using Linblad formulation 

M 

Furthermore, the thermodynamic properties of Davydov-Scott monomer con- 
tacting with thermal bath has also been investigated using Lindblad open quan- 
tum system formalism through path integral method. It has been found that the 
environmental effects contribute destructively to the specific heat, and large inter- 
action between amide-I and amide-site is not preferred for a stable Davydov-Scott 
monomer 

Unfortunately, although represented in non-relativistic hamiltonians, there is 
still no clue on how to expand further the hamiltonian. Because there is no basic 
assumption which is generally valid to restrict ourself when introducing the new 
terms of interaction. Therefore, we have proposed a new model based on the la- 
grangian with certain symmetry. The model was originally motivated by a need to 
take into account the effects from surrounding fluid, since the living systems are 
always immersed in a fluid medium. Such models are already available, one of them 
is the so-called fluid QCD representing the quark-gluon-plasma (QGP) system of 
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highly densed and strongly interacting plasma l^M The model has been ^apph ed to 



M 

describe the high energy QGP and plasma dominated early universe l lz l io [ The 
model is justified by the fact that it induces the non-relativistic Eule r equation 
with certain external forces when one takes the non-relativistic limit ^1. 

Concerning the fact that an (elementary) biomatter has no intrinsic degree of 
freedom like spin, it is considerable to represent its elementary constituents as scalar 
(boson) fields governed by the bosonic lagrangian, 

^matter - (^$) f + V(§) , (1) 

where V($) is the potential. For example in the typical $ 4 — theory, 

K($) = -im|$t$_^A(#$) 2 , (2) 

where m$ and A are the mass of matter and the dimensionless coupling con- 
stant of matter self-interaction. One can impose the above bosonic lagrangian 
to be gauge invariant under local (in general non-Abelian) gauge transformation, 
U = cxp[-iT a 9 a (x)} wl - iT a 9 a {x) with 6 a < 1. T a 's are generators belong to 
a particular Lie group and satisfy certain commutation relation [T a ,T b ] — if abc T c 
with f abc is the anti-symmetric structure constant. The matter field is then trans- 
formed as $ $' = cxp[— iT a 9 a (x)} $, with T a are n x n matrices while $ is an 
n X 1 multiplet containing n elements, i.e. 

and $ T = ($1 $ 2 • • • $„) , (3) 



4> 



$ 2 



for n dimension Lie groups as SU(n), 0(n + 1), etc. It is well-known that the sym- 
metry in Eq. ([T]) is revealed by introducing gauge fields A a ^ which are transformed 

as U£ U£' = U° - ^{d^6 a ) + f abc e b U<, and replacing the derivative with the 
covariant one, 2? M = 9 M + igT a U^. Anyway, the number of generators, and also 
gauge bosons, is determined by the dimension of group under consideration. For an 
SU(n) group one has n 2 — 1 generators and the index a runs over 1, 2, • • • , n — 1. 

The gauge boson is interpreted as a "fluid field" with velocity u^, and takes 
the form [7™ — (U§, U a ) = u®<fi with u M = 7 a (l, — v a ) and tj> is an auxiliary boson 

field, while 7° = (l — |v Q | 2 ) This is nothing else than rewriting a gauge 

field in terms of its polarization vector and wave function which represents the fluid 
distribution in a system. It has further been shown that the non-relativistic fluid 
equation can be reproduced 

For the simplest case of single matter to represent the single strand of DNA, 
using the Euler-Lagrange equation one obtains an EOM as follow ^21, 

(0 2 + m| + 2g 2 U 2 ) $ + ^A $ 3 = . (4) 
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for a real $ field. This can be solved analytically and the solution is well known. 
It can be concluded that the surrounding fluid enhances the amplitude at short 
distance ^31. Moreover, the double stranded DNA can be modeled using the same 
lagrangian but with n = 2. Again, the Euler-Lagrange equation in this case yields 
the following EOmQ31, 

(d 2 +m%- 4iga 2 U^d^) $ + ^yA ($ T $$) = , (5) 

at non-relativistic limit for constant fluid velocity and 4>. Here, 4> T = (fa, fa). This 
can be solved analytically under further assumption as done in H31. 

Anyway, so far we have concerned only the EOM's from the lagrangian. All 
of these would end up to the nonlinear differential equations which are in most 
cases unsolvable. Therefore one should consider another way to extract meaningful 
physical observables. 



3. Statistical mechanics observables 



For the sake of simplicity let us consider a simple example of protein dynamics using 
the so-called fa theory ^21 The calculation has been performed using statistical 
mechanics and path integral method. In particular, the evolution of heat capacity 
in term of temperature has been investigated for various levels of the nonlincarity 
of source, and the strength of interaction between protein backbone and nonlinear 
source. It has been argued that the nonlinear source contributes constructively to 
the specific heat especially at higher temperature when it is weakly interacting 
with the protein backbone. This may indicate increasing energy absorption as the 
intensity of nonlinear sources are getting greater. 

Similar to the pr evio us case, the model adopts the fa lagrangian but with ad- 
ditional interactions ^1 



Ctot = £ C (<P) + C s (ip) + c int (4>, i/j) , 



where, 



C = r 



(d,fa f (d»<t>) - 



A (0t ) tyty) 



(6) 

(7) 

(8) 
(9) 



representing the conformational changes of a protein backbone and the nonlinear 
source injected to the backbone, while the last one is the interaction term between 
both. Imposing a local U(l) symmetry to the total lagrangian and considering its 
minima lead to the vacuum expectation value (VEV), (ip) = y2A/A (</>). This non- 
zero VEV then yields the so-called spontaneous symmetry breaking. The lagrangian 
yields two EOM's in term of 4> a- n d ip respectively. However we are not interested in 
discussing its dynamical behaviors. The statistical observables can be conveniently 
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calculated from the generating functional by the perturbation method The gen- 
erating functional for scalar fields is written as, 



Z = I V(j>T>ip exp 




drdx Ctot (0, VO 



(10) 



at the finite temperature case in Euclidean coordinates. The integral in Eq. (fT0| can 
be evaluated analytically using the Gaussian integral. This can be accomplished by 
rewriting it in term of Gaussian integral using the Fourier representation of Green's 
function to obtain the master equation 



iVexp 



i Jo Jo 



dx 



-~AA$(0)+AA*(0)A*(0) 



(11) 



For statistical observables, Let us consider the specific heat of the system in 
a constant volume, Cy, that is a particular interest from experimental point of 
view. The specific heat can be derived directly from the partition function using the 
relation, 



c v = fr 



d 2 hxZ 

dp 2 



In the present case, it has been found to be 



EE 



d 2 



-aa2,(o)-aa^(o)a^(o) 



(12) 



(13) 



after performing the integration over r and x respectively, while the overall factor 
N = [47rsinh(fc/3/2)] _1 . 

The integration in Eq. (|13[) should be done numerically. It has been argued that 
both A and A contribute to the evolution of Cy in an opposite way, and could 
completely cancel each other at certain values. This occurs when the symmetry is 
maximally broken which means that increasing energy absorption prefers high level 
of nonlinearity of sources and at the same time weak interaction between the sources 
and protein backbone. 



4. Conclusion 

New models for biomatter dynamics using lagrangian have been presented. The 
approach have been applied to describe typical cases of single and double stranded 
DNA, and the protein folding mechanism. The usage of statistical mechanics method 
to extract the specific heat has also briefly been demonstrated. It can be argued 
that the approach might be able to provide new directions on modeling the living 
phenomena at elementary level. 

However more theoretical works are still going on. More importantly it should 
be justified by experiments in the future. 
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